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Abstract 

, The Stokes equation with the varying viscosity is considered in a thin tube structure, i.e. in a con- 

nected union of thin rectangles with heights of order e << 1 and with bases of order 1 with smoothened 
boundary. An asymptotic expansion of the solution is constructed: it contains some Poiseuille type 
flows in the channels (rectangles) with some boundary layers correctors in the neighborhoods of the 
bifurcations of the channels. The estimates for the difference of the exact solution and its asymptotic 
approximation are proved. 
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(N 

(N ; 1 Introduction 

i— i ■ 

The blood circulation, the transport of cells and substances in the human body as well as some liquid- 
. cooling systems and oil-recovery/oil-transport processes in engineering, are modeled by the equations of 

fluid motion posed in thin domains. The present paper studies the Stokes equation with varying viscosity 
in a tube structure. In two-dimensional case a tube structure^] (or pipe-wise structure) is some connected 
union of thin rectangles with heights of order e and with bases of order 1 with a smoothened boundary 
([5]. jlOj): here e is a small positive parameter. An asymptotic expansion for the case of a constant 
viscosity have been constructed in [1] . It is " compiled" of some Poiseuille flows inside the rectangles glued 
by some boundary layer solutions, in the neighborhood of the junctions. Here we consider a more general 
case when the viscosity is not constant but depends on a longitudinal variable for each rectangle. This 
situation models a blood flow in a vessel structure where the viscosity depends on the concentration of 
some substances diluted in blood or some blood cells. Indeed, the asymptotic analysis of the convection- 
diffusion equation set in such domains ([3], |10j ) shows that in the case of the Neumann (impermeability) 
condition at the lateral boundary and small Reynolds numbers, the concentration is asymptotically close 
to the one-dimensional description, that is the convection-diffusion equation set on the graph. The solution 
of the problem on the graph is the leading term of the asymptotic expansion, and it evidently on depends 
on the longitudinal variable. On the other hand, the viscosity often depends on the concentration of the 
diluted substances or distributed cells, and so, it depends on the longitudinal variable. Of course, the fluid 
motion equation is coupled with the diffusion-convection equation in this case. However, if the velocity 
is small (in our case, it is of order e ), then neglecting the convection, in comparison with the diffusion 

1 see Fig[5] 
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term or iterating with respect to the small term @, we get the steady state diffusion equation; in absence 
of the source term in the right hand side, it has a piecewise-linear asymptotic solution on the graph for 
the concentration. So, in this simplified situation, the diffusion equation can be solved before the fluid 
motion equation, and we obtain for the flow, the Stokes or Navier-Stokes equation with a variable viscosity 
depending (via concentration) on the longitudinal variable. These arguments provide the motivation for 
the Stokes equation with variable viscosity; to our knowledge, it has not been studied earlier from the 
asymptotic point of view. There are, of course, many other practical problems involving fluids with variable 
viscosity. For example, the presence of bacteria in suspension (see [3]) may change locally the viscosity. 

In the first part, we consider the case of a flow in one rectangular channel with the periodicity condition 
at the end of the channel. An asymptotic expansion of solution is constructed and justified (a regular 
ansatz). 

In the second part, we consider the case of a flow in one rectangular channel with the inflow/outflow 
boundary conditions at the ends. In this case as well we construct an asymptotic expansion of solution, 
that contains a regular ansatz and the boundary layer correctors. 

In the third part, we construct an asymptotic expansion for solution of the Stokes equation set in a 
tube structure. Here as well we associate to each rectangle a regular ansatz and then we glue all the 
regular ansatzes with help of the boundary layer correctors exponential decaying from the junctions of 
rectangles. This procedure is similar to the procedure described in [lOj . 

In the fourth part we consider some numerical experiments comparing the numerical and asymptotic 
solutions. 

All constructed expansions are justified by calculation of residual terms and application of a priori 
estimates. 

2 Flow in one channel 

Consider a small parameter e, e = |, q G N, and define then a thin domain 

D £ = |(xi,.t 2 ) G K 2 : < xi < 1, -| < x 2 < || • 

Assume that incompressible, viscous fluid fills the domain D £ . Let / be the exterior force applied to the 
fluid. 



Figure 1: Thin domain 



2 The simplest coupled diffusion-convection equation is 

-div(u(C e )Vu £ ) + Vp e = f(x) 
divu e = 0, 

-AC e +%.VC E = 
where C E is the concentration, u e is the fluid velocity and p e is the pressure. 

If u s is small with respect to C e , we can write an expansion for C £ with respect to the small parameter that is the ratio of 
magnitudes of |u e | and |C e |. Then for the terms of this expansion, the third equation may be solved before the fluid motion 
equations. Another possible approach is the successive approximations (fixed point iterations) : 

-div(Kci n) )X>4 n) ) + Vpi n) = f{x) 
div4 n) = 0, 

-Aci n) +4 n-1) .Vci n) =0 

where n is the number of the iteration. 

In both approaches, we get a problem with the variable viscosity. 
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Consider the following steady state Stokes problem: 



(1) 



' -div(p(x 1 )Vu £ )+Vp e = f(x) inD e , 

divu £ = in I> £ , 

u e (x\, f ) = for x\ e (0, 1), 

u £ (xi,-f) = for xi 6 (0, 1), 

u £ {Q,x 2 ) =eVo(f ) for s 3 e (-§,§), 

k u e (l,a:2) =eVi(^) for x 3 e (-§,§), 



The unknowns of this system are the velocity u e and the pressure p £ of the fluid. 

The non homogeneous boundary conditions for the velocity are given by the functions ipo , ipi whose 
second component is equal to zero, and the first satisfies 

i i 

(2) j V?oi(&)d6= /%ii(6)d£ 2! 

where G Co([ — |, |]). Here (2>it e )y = \ + i^f 2 "^ v satisfies the following conditions 

(3) v(xi) = v + vi{x\) 

where v\ g Cq([0, 1]). Moreover there exist p > 0, k > such that v(x\) > n for all x\ G (0,1) and 
^i(xi) = for all xi G (0, p) U (1 — p, 1) (and so, > 0). 

2.1 Variational formulation of the problem 

In order to obtain the variational formulation of problem (TTJ), let us introduce the following space 

(4) H(D £ ) = {Hl(D e )) 2 : div(u) = 0} 

and assume that / G (L 2 (D £ )) 2 . 

Applying the extension theorem (see [5]), one can find a function ip e G (if 1 (I? e )) 2 such that divip £ = 0, 
ipf = ipo and ipf = if i and ipf = . Changing the unknown function u £ by v £ = u £ — e 2 ip £ , 

we give the variational formulation of problem (TTJ): 



(5) 



/ v(x\)Vv e : Vip = f f-^-e 2 \ v{x x )Vip e : Vip, \/ip G H(D £ ). 

J D e JD E JD e 



Definition 2.1. We say that u £ is a weak solution of problem (TjTJ) if v E — u £ — £ 2 </r G H(D £ ) and satisfies 
©• 

Proposition 2.1. If u £ is a weak solution for problem |T]) ; then there exists a distribution p £ £D (D £ ) 
such that (u e ,p s ) satisfies (TT])i this problem in sense of distributions. 

Proof. If we take v £ = u £ — e 2 ip £ , then from (5) we have that 

(div(i/(xi)©w e + e 2 ^(xi)P^ £ ) + /, ^) = 0, W e (D e ) 

From De Rham lemma, it follows that there exists a distribution p e , unique up to an additive constant, 
such that 

-div(i/(xi)X>?; e ) + £ 2 d\v{v( Xl )Vp £ ) - / = -Vp. 
Theorem 2.1. The variational problem ([5]) admits a unique solution v £ G H{D £ ). 

Proof. The Riesz's theorem gives the existence and the uniqueness of solution because the norms 
||u £ ||j = J J v(x 1 )Vv £ : Vv £ and \\v £ \\^ h i( De)) 2 arc equivalent. 

As a consequence, we have: 

Proposition 2.2. Let v £ G H(D £ ) be the solution of the variational problem Then the following 

inequality holds 

\\v e \\(h(d s ))* < C(n,C £ PF ) ^||/|| (i,2(n e ))2 + \\¥> £ \\(Hi(D e )) 2 ^ > 

where Cp F stands for the Poincare-Friedrichs inequality constant and K is the lower bound of the viscosity 
©■ 
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Mention that Cp F can be estimated by eC, where C is independent of e (see [TU]). 
Consider the case when (p e = 0. Then we have: 

Proposition 2.3. The following inequality holds 

\\VPe\\ H -HD e )<C\\f\\ (L2{De)r , 

where C is a constant independent of e. 

Proof. From De Rham lemma, it follows that there exists a distribution p £ , unique up to an additive 
constant, such that 

V Pe = div(i/(a;i)XHt e ) + /. 



It means that 



/ p £ div^= / vfayDue-.Thl)- I /•</>, WeH(D £ )) 2 , 
Jd, Jd, Jd, 



and so, 

\S D Pediwip\ \J D v(x 1 )'Du e : Thp - J D f ■ 
sup j. j. = sup < 

< (C(K,Cpp) + Cp F )\\f\\( L 2( De )y, 

where C(k, C pf ) is the constant of the a priori estimate for u £ (see above) and Cp F is the Poincare- 
Friedrichs constant (it is of order e). This inequality gives the estimate of Vp e in the if -1 — norm. 

An asymptotic analysis of problem ([1} shows that an asymptotic solution is given by a Poiseuille 
type flow, with two boundary layer correctors localized in some neighborhoods of the ends of the channel. 
Didactically, it would be better to separate the construction of the Poiseuille type flow for varying viscosity 
and the construction of the boundary layer correctors. That's why we simplify problem ([T]), and replace 
©55 ©6 by the periodicity condition with respect to x\. 

Introduce the Sobolev space 



H per (D e ) = {ue (H^ erth0 {D e )) 2 :divu = o} 



Here Hp er 1 {D e ) is the completion (with respect to the 7J 1 (D e )-norm) of the space of the C°°(Rx [— |, |])- 
functions vanishing at the boundary x% = ±f and 1-periodic in x\. As in the beginning of the section, 

/ e {l\d £ ))\ 

Definition 2.2. We say that u £ £ H per (D E ) is a weak solution of the periodic problem 

—c\\\{v(xi)T>u e ) + \7p £ = f(x), in D £ , 
divu £ =0, in D £l 

(6) < u e {xi,l) =0, forx x e (0,1), 

u e (xi,-§) = 0, for xi 6 (0,1), 

u £ is 1 — periodic in X\, 

if and only if it satisfies the integral identity 

(7) f v(x x )Vu e :V^= [ f.t/i, G H per (D s ). 

JD C JD e 

As in theorem 12. 1[ we apply the Riesz theorem and prove the existence and the uniqueness of u £ G 
H per (D £ ), solution of problem The a priori estimate is given by 

\\ue\\(H l (D e )y < C(k, C P p)\\f\\ {L 2 {De)) 2 

where C(k, Cp F ) = O (e) and k is the lower bound of the viscosity ([3]). 

Proposition 2.4. If u £ is a weak solution for problem © then there exists a distribution p £ € T> (D £ ) 
such that (u £ ,p £ ) satisfies (|6]) i in sense of distributions, and 

l|Vp e || B -i (D .)<C||/|| (L3(D . ))a , 
where C is a constant independent of e. 
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Figure 2: Infinite layer 



Proof. By J7J 



v{x\)Du e : T> ip 



f-iP, Vi'€H per (D e ) 



From De Rham lemma, it follows that there exists a distribution p e , unique up to an additive constant, 
such that 

Vp e = d.iv{v(xi)Vu e ) + f. 

It means that 



p s dbri> = / u{ Xl )Vu s :V^- f ■ $ (i/„ 1 er ,i,o(^)) 2 , 



and so, 



sup 



I In PedivV'l 



sup , , ,, < 



^(H^ r ,i,o(D c )V ll V ^ll(L2(D e ))4 V'e(^ rjl ,o(^)) 2 IIW||(L*(A0)* 

< (C(K,Cpp) + Cp F )||/||(_ £ ,2( £ ) e ))2, 

where C(n,C PF ) is the constant of the a priori estimate for u e (see above) and Cp F is the Poincarc- 
Fricdrichs constant (it is of order e). This inequality gives the estimate of Vp £ in the H~ x — norm. 

2.2 Asymptotic analysis of the problem 

Let us first construct the asymptotic expansion for the solution of the periodic problem ([5]); then we will 
study the non periodic one. 
Define the infinite layer 



{ (2:1,2:2) 



x\ e R, -- < x 2 < 



1} 



Assume that / = /iei, f\ 6 C°°(R), fi is 1-periodic in x\. 

Denote {ip) 1 = J Q ip(xi, ai2)da;i, (-0) 2 = Xfi ip(xi,x 2 )dx2- An asymptotic solution is written as: 



(8) 



u\{xi 1 x 2 ) 
u\{x x ,x 2 ) 
p k {x ll x 2 ) 



3=0 

k 



X2 

e 



J2 £J 



3=0 

Si 



3=0 



3=0 



with ui^,U2.j,Pj and gj 1-pcriodic functions in 2:1, such that (pj) 2 = 0. 

Substituting ([5]), in ([5]), equating the coefficients of the same powers of e and denoting ( 2 = 7 we 
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obtain: 



f -afrOi)%f 



(9) 



t/(xi) d 2 Ui j _ t/(xi) d 2 u 2 ,j-2 _|_ dpj-i , dq 3 _ r r 
2 2 a^O^i 9*1 + 3n — /l°jO) 



1 a 

2 9xi 



9^1 j , 9u2 



0. 



9a; i 3£ 2 

til,j (Xl,±i) = 

, U2, 3 (xi, ±|) = 



Denote: 



(it satisfies A^' = 1, -/Vi(±|) = 0) and 



^2(6 



2 I 42 4 



6 



iVi(r)dr, 



(here AT 2 (±; 



Aq(-r)dr = -^). Denote: 



D- 1 : F 
D- 1 : F 
D- 2 : F 



6 



F(x 1 ,r)dT, 



6 



F(xi,r)dr- 



F(xi,r)drd0, 



€2 /■« /i 
F(x X) r)dTdfl- (6+0 



F(xi,T)drd6». 



Theorem 2.2. The unknowns of (J^J) Uij,U2j 1 Pj,qj are given by the following relations: 



(10) 



Ml 



U2, 



-D 



|^^+ 2 f « f^( xt)^-^ - 

•° 1Z) 2 { a^a^f + S^T wtl7 (ail ( v ( 
r (ft ~ W) JV*(6), 



xi ) 9ui - j - 

■ Ll > dx! 



dp 
dx 



Proof. Integrating twice and using boundary conditions (jl]) 4! we get ([lOp,. This relation gives 
an expression of the unknown u\ .j via qj. All other functions contained by this relation are either known 
from previous computations or equal to zero. We integrate next the incomprcssibility condition §§§ 3 with 
respect to £2 with the boundary condition® 5 and get U2 j . Finally, integrating ^ 2 we get pj. The 
unknown function qj is determined from the boundary condition U2j (xi, ^) = 0. Actually, the boundary 
conditions ([9]) 4 and U2,j (xi, — \) — arc satisfied by the definition of D~ 2 and D , while U2.j (xi, 5) = 
gives: 



(11) 







1 



- M 



30 



dx\ \6i^(xi) \dx\ 



D~ D~ 



d ( .Buxa-i 



\ d&dx 2 dxi \i^(xi) \dxi \ y " v '" i/ Q Xl J dx\ 
In particular, for j ~ 0, we have the Darcy equation for the leading term of the pressure qo~. 

< 12 > ST (sfcj (&-*))-<>■ 



G 



and so go is the 1-periodic function given by the formula 



,.<*,) = f (aw - f^M) dx - (J" (aw - f^w) 

Then ui, (xi,&) = ^ (fir ~ /i) ^1(6), «2,o(zi,6) = and po(<ci,60 = 0. 
For j = 1 q 1 is the 1-periodic solution of equation 



d f 1 ggi 

It follows: qi = 0, «i, 1(^1,^2) = and ^2,1(^1,^2) = 0. For j = 1, (fTUjl o yields: 



For j = 2, (|TTj) becomes 



(13) 75- 



1 <9g 2 



(9 



3 / <9<?i 



9xi V9xi 



= 0. 



-h (iVi(6)-(iVi) 2 ). 



\6^(xi) 9xi 

where D~ 1 D~ 2 (-/Vi(£ 2 ) — (-^1)2)1 = 7\L 7^ 0, g 2 is 1-periodic solution of given by: 



«2 = 



«2 = * 



- 1 



ft(a?i) = -12^(/i)i 

It follows: 

< 

t «2,2(^1,6) = (f^ - /l)] ^D" 2 (^(f ) - (iV 1 ) 2 ) - ^ (^^(f)^ 

and we get the following theorem. 

Theorem 2.3. The asymptotic solution of problem is given by: 



(14) 



u§(xi,x a ) = -e 5 (3§7 



2 a 2 



p fe (xi,x 2 ) 



'iu 



afr"(Klr)^(t)lft))+ ( £6 )' 

+ ^ 2 + sir (fir - A) ) - + °( £3 )' 



where qo and q 2 are periodic solutions of obtained as the unique solution of (jlip (| 13[) respectively. 
Remark 2.1. Solution given by (|14j) can 6e written as follows: 



u\{x llX2 ) = -^^^N 1 (^)-eU^^u''(x 1 )D- 2 (N 1 (f)-(N 1 ) 2 ) 



2 </i)i. 



24TC(/i)i / i/ (xi} \/at 1 , /)f,6\ 

+ (Kxo)! (ti^tt) iV 2(-) J + <a £ 



- f ( AW - ir" M ) dT - (f ( /1(T) - fr" <T) ) dT ), 

Mi 



(-12^(/i 



^1/ ( 2;i )(iVi(f)-(iV 1 ) 2 )j+0( £ 3 ) 
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Introduce the following function: 



F k {x u x 2 ) 



2- ( vfxx)^^ ) + v ^> " "'u*- 1 - SEt ) ei 

Xi V V ' 9a!! y 2 02:1042 crai y 1 



k(xi) d 2 u 2 , k - 



(15) 



1 a 

2 9xi 



'(*0 

Theorem 2.4. Let (u k ,p k ) be the asymptotic solution given by © and (u,p) the solution of ©• Then 
the following estimate holds: 

HW = 0(e^i), 



\u k -u\ 



Proof. Denote (U ,P ) = (u — u,p — p). We obtain the following problems for (U ,P ) 



(16) 



' -div(i/(xi)W c ) + VP fe = -£ k + 1 F k (xx,x 2 ) in D e 
divU k =0 in D e 

U k {x u %) = for Xl G (0,1), 

C/ fc (a; l! -f) = 

[7 fe is 1 — periodic in x\. 



for xie(0,l), 



The a priori estimates give us: 



u k -u\\ m{De? <C{K,C PF )\\e k+ "F k \\ LHDc) ,=0[e 



M+ 1 . 



and 



\Vp k -Vp\\ H - 1{De) = 0(e 



Let us consider now the non periodic case, i.e. let us construct an asymptotic expansion of the solution 
for problem ([1]). Assume 

f = h{x 1 )e 1 , /iGC°°([0,l]). 

Define an asymptotic solution by: 



(17) 



u k (x u x 2 ) = u k (x u x 2 ) + u k BL0 ^)+u k BL1 (^,f ) 
p k (x u x 2 ) = /(x 1 ,x 2 )+p| L0 (f)+p| il (^,f ) 



The expressions of u k , p k are the same as above: flSJ) , (fTU)) , (fTT1) 

Since the functions given by ([5J do not satisfy exactly the boundary conditions ([I]) 5 , (H]) 6 we add the 
boundary layer correctors. They correspond to the left end for i = and to the right end for i = 1 and 
their expressions are given by: 



(18) 



k 

3=0 



Pbl 



X\ - l X2 

e e 

Xi — i X2 

£ £ 



Here i = 0,1. To obtain problems for the boundary layers corresponding to the left side, we introduce the 
domain 11+ = (0, 00) x (— |, The problem 



(19) 



f -$A^(0 + V tf §(0 = 
div£u° = 



,,0/n t \ - ln X ( u l,j(°>60 

. ^ (0 ' 6) -^ ~U^-i(o^ 2 ; 



ifCen+, 
if£en+, 

if6 6(-i,i), 



with the compatibility condition 
(20) 



(VOl^'O - Wl,i(0,6)} 2 = °- 
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where (■0), = /Ji ^(^j £ 2 )d£2, will give the boundary layer correctors for the velocity and for the pressure 

2 

corresponding to the left end (see [4]). This condition (|20| generates a boundary condition for qf. 



(Voi) 2 <fjO + 



6z/ n 



9gj(0) 
dxi 



-Wi(o) + 



_ 2 d 2 U 2 J-2 



2ir 



9x 2 



^0 9a; 1 



= 0. 



1*1=0 / 



In a similar way we introduce the boundary layer correctors corresponding to the right side. The boundary 
layers for the velocity and pressure are defined on II = (— 00, 0) x (— |, |). The analogous boundary 
condition for qj is satisfied automatically because of the conservation of the average of u\j (xi , £2) ■ Indeed, 

(«i,j(i»i,6)>2 = ^"V,.? (31,6)1, _i ; 

on the other hand, the condition 

«W| ta . 4 =-^- l «i J (*i,6)| <a . i =o 

is equivalent to the equation on Oj, i.e. this equation on Oj is equivalent to the conservation law for the 
average {u± t j(xi, £ 2 )} 2 - The problem satisfied by the asymptotic solution of order k is as follows: 



(21) 



-div(v( Xl )Vu k ) + Vp k = f -e k+1 F k (x u x 2 ) 
divu k = 

U e (xi, |) = 

u £ (xi, -§) = 

it 

u e (0, x 2 ) = e 2 Mf ) + £ k+ Wk(®, 6)e 2 + £V 1 

j=o 
k 

fi e (l, ia) = eVi( f ) + £ fc+3 u 2lfc (l, 6)e 2 + XV +2 u§ 



J+2J1) 

j 



1 x 2 



(0) / 1 x 2 



in £> e , 
in D e , 

for xi G (0, 1), 
for xi G (0, 1), 

for x 2 G (-§, |), 
for x 2 G (-§, |), 



Here F fc is given by (TT51) . 

Mention that the boundary conditions for u fc on x\ =0,1 are not yet satisfied exactly because the 
traces of each boundary layer on the opposite side of the rectangle are exponentially small but do not 
vanish completely. These traces should be eliminated by a small additional corrector. 

Let us construct a new function U k which satisfies the same boundary conditions as u k on xi = 0, 1. 
Let us describe its construction: let U k : D e — > R 2 be a solution of the following problem: 



(22) 



' U k e{HHD £ ))\ 

divU k = 
U k (xi,l) = 
C/ fe (x 1) -|) = 

k 

U k {0, x 2 ) = £ fc+ Vfc(0, 6)ea + ]T 

3=0 

U k (l 7 x 2 )^e k +h l2 . k (l,^)e 2 +J2 £j+2u f ) (->- 

3=0 



in D £ , 

for xi G (0, 1) 
for xi G (0, 1) 



1 X2 

£ ! £ 



for X 2 G (-§, §), 

for x 2 G (-§, |). 



Proposition 2.5. Problem (|22[) /ias ai Zeasi one solution with the property 



(23) 



ll C/fc |l(ff 1 (£ , e)) 2 



Of£ fc+ ^ 



Proof. Define : D 



I 2 , where D = (0,1) x (-|,i), 



(2/1.2/2) 



e{U k ) 
(U k ). 



with (yi, 2/2) = (a-i, — ). 

£ 
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Obvious computations lead to the following problem for w k 



div^to* = 
w k { yi ,-\)=Q 



w k (Q, y 2 ) = e fc+3 M 2 ,fc(0, y 2 )e 2 + 



(24) 



( (-In 



3=0 
k 

E 

3=0 
k 



• >)2 



in D, 

for yi G (0,1), 
for yi G (0,1), 



for y 2 G (-i, |), 



= e fc+3 M2,fc(l,y 2 )e 2 + 



E^S 7 



k 

E 

\ 3=0 



,2/2 



for y 2 G (-|,|), 



/ 



Here div y is the divergence in variables. As in [5] we can prove that there exists a function w k G 
( y H 1 (D £ )) 2 so that 1 1 to* || (jyi(zj)) 2 < C\\w k '\\ ✓ ^ (g^ 2 wrtn ^ e constarL t C independent of e. Using the 
properties of the boundary layer correctors (their exponential decay rate), we get: 

IK fe ||(HW = 0(e k+s ) 

Direct computations give ||t^ fc ||(#i(£> e )) 2 < ~xll' u; e \\(h 1 (d)) 2 - Combining these two estimates we achieve 
the proof. 

The function 



(25) 



U k = u k - U k 



satisfies the same boundary conditions as it in x = 0, 1. The problem for the new functions U k , p k is an 
obvious consequence of (|2"2"j) and ([2"Tj): 



(26) 



-div(v( Xl )VU k ) + Vp k = f - e k+1 F k ( Xl ,x 2 ) - &\v{v{ Xl )VU k ) in D £ , 

divU k = in D e , 

&(*!,§)= for Xl G (0,1), 

U s {x u -f) = for Xl G (0,1), 

U e (0,x 2 ) = eVo( 2 r) for ar 2 G (-§, §), 

[ C/ £ (l ) x 2 )=eVi(t) for x 2 G (-§,§). 



Theorem 2.5. Let (u k ,p k ) be the asymptotic solution given by (fl~7|) and (u,p) the exact solution of (JTJ) . 
Then the following estimates hold: 



(27) 



I" - u e\\H(D e )- 

\Vp k -Vp £ \\ H - 



= e k + 



= O 



(>+§) 



Proof. From (25} it follows: 



^|| ff (D £ ) 2 < 11^" - u e\\ H (D E 



\W k \\ H 



(De) 2 



= 0[e k+ i 



0{e k+ i 



From (fT5" ]), (j2"5|) and (|2"6"|) it follows: ||u fe — ull^-^^ = O ye k+ ^j. The estimate for pressure is a consequence 
of ([2T| 1 and the a priori estimate. 



3 Flow in tube structures 

In this section we are going to construct an asymptotic expansion to the solution of problem (JT|), stated 
in a tube structure containing one bundle. We justify the error estimate. Let us define a tube structure 
containing one bundle. 
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Figure 3: One bundle of segments B 




Figure 4: The rectangles Bj 

Let ex, &%, ■ ■ ■ , e n be n closed segments in M 2 , which have a single common point O (i.e. the origin of 
the co-ordinate system), and let it be the common end point of all these segments. Let fa, fa, ■ ■ ■ , fa be 
n bounded segments in M 2 containing the point O, the middle point of all segments, and such that fa is 
orthogonal to ej (for simplicity assume that the length \fa \ of each fa is equal to 1). Let /3| be the image 
of fa obtained by a homothetic contraction in ^ times with the center O. Denote Bj the open rectangles 
with the bases fa and with the heights ej , denote also faj the second base - side of each rectangle Bj and 

let Oj be the end of the segment ej which belongs to the base fa- (see Fig. [4]). Define the graph of the 
tube structure as the bundle of segments ej having a common point O (see Fig. [3]) 

n 

B=\je j . 

3 = 1 

Denote below Oo = O. Let 7?, j — 0,1,..., n, be the images of the bounded domains jj (such that jj 
contains the end of the segment Oj and is independent of e) obtained by a homothetic contraction in | 
times with the center Oj. 

Define the tube structure associated with the bundle B as a bounded domain (see Fig. [5]): 

-((usf)u(lH ? ))' 

Here the prime stands for the set of the interior points. Assume that dB £ G C 2 (the result may be 
generalized for the case of the piecewise smooth boundary dB £ with no reentrant corners). Assume that 
the bases f3 £ of B e , j = 1, ...,n, arc some parts of dB e . We add the domains 7?, j = 0,1,..., n, to 
smoothen the boundary of the tube structure. 
Consider the following system of equations: 

( -&iv(v(x)Vu e ) + Vp E = f(x) in B £ , 
(28) < divu E = in B e , 

y u e = g on dB e . 

Here, g = on the lateral boundary of the rectangles composing B e \ moreover 5 = anywhere with the 
exception of the sides /3| of the rectangles Bj (these sides are assumed to belong to the boundary of the 

tube structure); g £ C 2 (/3|), and for each j, g — £ 2 gj( x ~ e 3 ) on faj, the vector valued functions gj £ C 2 
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Figure 5: One bundle tubular structure B e 



do not depend on e. Let / be a vector-valued function of (L 2 (B £ )) 2 . The solvability condition gives the 
relation 



(29) / g.nds = 0. 

Introduce the local system of coordinates Ox\ 3 x e 2 3 associated with the segment tj such that the direction 
of the axis Ox\ 3 coincides with the direction of the segment OOj, i.e. x\ 3 is the longitudinal coordinate. 
The axes Ox\ 3 x 2 3 form a Cartesian coordinate system. Denote doe the infimum of radius of all circles with 
the center O such that every point of it belongs only to not more than one of the rectangles B e , j = 1, . . . , n 
and di is the maximal diameter of the domains 70, 71, . . . , 7 n . We finally introduce the notation 

doe = max{(io£, die}. 

Consider the right hand side vector valued function / " concentrated" in some neighborhoods of the nodes 
Oj and diffused in the rectangles, i.e. 

/ = $j ^ - £ ° 3 ^j , {oi\x-O j \<d e,j = 0,...,n, 

f = for \ x - °j\ > d o£, x\ j G (0, \ej\),j = l,...,n. 

Here f 3 G Co°([0, |ej-|]),$j G C^(Q), (j = 0, 1, ...,n), where Q is a ball |£| < d . Assume that 

(31) v{x) = vo + Vj{x\ 3 ) 

such that Vj{x\ 3 ) = for all x\ 3 G [0,/3] U [|ej| — /3;|ej|], where (3 is a positive constant such that 
P < mirij^-; v G C 2 and there exist ko G M + such that v(x) > for all x G £> e . Without loss of 
generality we may assume that /j (2^) = for all x\ 3 G [0,/3] U [ | e ^ | — /3; |ej|]. 

Let Hdi V= o{B € ) be space of the divergence free vector valued functions from H 1 (B e ). Let H^ iv=0 (B e ) 
be the subspace of vector valued functions of -ffdiv=o(S e ) vanishing at the boundary. Assume that g can 
be continued in B e as a vector valued g of Hai V= o(B £ ) . The variational formulation for ([2"5f is as follows: 
find u e G Hdiv=o(<6 e ) such that v e =!i E -j£ H^ iv=Q (B £ ) , and such that it satisfies to the integral identity 

/ y(x)Vv e :V<p= [ f ■ tp - [ v(x)Vg : Dtp, G H° iv=0 (B £ ). 

JB E JB' JB e 

The Riesz theorem will give the existence and the uniqueness of such a solution because the norms 
||u|| = J j Bc v(x)T>v : T>v and || f|| (fyi (B«)) 2 a re equivalent. We have as a consequence that 

II^Hcff 1 ^)) 2 < C(C PF , Ko) (||/||l 2 (£^) 2 + 11.911 (/J 1 , 

where Cpf is independent on e (see [10]) and kq is the lower bound of the viscosity (I3T1) .. 
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Proposition 3.1. If u £ is a weak solution for problem (|28[) then there exists a distribution p £ £D (B e ) 
such that (u e ,p e ) satisfies (|28p i in the sense of distributions. The following inequality holds in the case of 
9 = 0: 

||VjD e || ff -i(B=) < C||/|| (L 2 (Be)) 2. 

Here C is a constant independent of e. 

The proof is similar to that of the Propositions of section [5J 

3.1 Asymptotic expansion 

We construct the main part of the asymptotic expansion in a form 
(32) M " = ]TV+ 2 j J2 u^ L )x,{x) + Y J ^( X -^ 

1—0 I e— ej l,...,n 2—0 



(33) 



1=0 I e=ej;j=l,...,n i=0 ^ ' J 

k n 

e—ej-.j — l : ...,nl—Q i— 



IS 



.T e ' L = (a;^, ^-). Later, in the end of the section we will add an exponentially small corrector multiplying the 
boundary layers by a cut-off function 77 in the subdomain where the boundary layers are just exponentially 
small (see (|42l) . (T4"3"|) ). The last sum in (|2"5)l is taken for all nodes Oi and the value (x e { = 0) is calculated 
at the point x = Oi, function q^ is supposed to be continuous on the graph B. Here Xe( x ) is a function 
equal to zero at the distance less than (do + l)e from Oj, j = 0,1, ... , n, equal to zero on the rectangle 
Bj if x\ 3 < (d + l)e or if — \ej\\ < (do + l)e; we suppose that function Xe is equal to one on this 

rectangle if x\ 3 > (do + 2)e and \x\ 3 — \e$\\ > (do + 2)e, and we define Xe by the relations Xs( x ) = X 

if (do + l)e < x\ 3 < (do + 2)e and ^(s) = x ( X ° J ; M ) if (d + l)e < |e,-| - a:? < (cf + 2)e. Here x i 
a diffcrcntiablc on K function of one variable, it is independent of e, it is equal to zero on the segment 
[—(do + 1); (do + 1)] and it is equal to one on the union of intervals (—00, —(do + 2)) U ((do + 2), +00). 
Moreover Xe is equal to zero on every 7J. The functions Oi, i = 1, n, are defined as follows 

leA 

9j = \x — Oj \ > mini—^- 
9j = l \x - 3 \ < irnj-^- 
The relation between the vector-columns x T and x ej,T (here T is the transposition symbol) is given 

by 

X T = TjX e i' T + O j = l,...,n 

where Tj is an orthogonal matrix of passage from the canonic base to the local one. Then applying the 
results of section [2J for every channel S|, we get uf, pf and qf defined up to the scalar constants cf, 

dq 3 

df. Indeed, denote q, 3 (x\ 3 ) the solution of equation (fTTj) with e . (0) = 0. Then the general solution of 

0x1 

equation (TTTj) has a form: 



(34) q? (x\ 3 ) = q? (x\ 3 ) + c\ 3 / v(s)As + d\ 3 , 

Ja 

where c\ 3 and d^ 3 are the undetermined constants; 
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where the second component of does not depend on c^ 3 , d^ 3 (see tUIlb); the same property holds for 
p\ 3 (see (Unjh); the first component (see (fTU]) i ) depends on c\ 3 : 



~e 
U 



3 — — D~ 
id ~ 



2,1-2 







v{x\ 3 ) 



^1,1-2 

dx\ 3 



hS j0 \+2c1 3 N 1 {C 2 3 ) 



d p?-i 

dx\ 3 



For x\ 3 £ [0, /?] lj[|ej| — /?, |ej|], , p^ 3 and dcfr 3 /dx\ 3 may be taken equal to zero because the right hand 
side fj is equal to zero for this values of x\ 3 , while the flow rate fp u^ 3 ^^ 3 )^2 3 i s constant on ej. Ni is 
the function introduced in section 2. 

To get the problems for the boundary layers, we introduce the domain Qo = U™ =1 f2j U 70, where fij 
are the half-infinite strips obtained from Bj by infinite extension behind the base pj and by homothctic 
dilatation in - times (with respect to the point O); let fij be obtained from Clj by a symmetric reflection 
relatively to the line containing /3j and let fioj = &j U 7] , where 7* is obtained from 7^ by a translation 
(such that the point Oj becomes O). 

Since v(x) = vq for all x\ 3 £ [0, [3] U [ | e ^ | — (3; \ej\], the boundary layer solution is a pair constituted of 
a vector valued function uf LOj and a scalar function pf L ° 3 satisfying to the Stokes system: 



(35) 



$o(0^,o- 



{cf (^oA ? (x,-^ (A^ 2 e ), 0) T ) - V« (Xi(^)^i)) - df +1 v a ^t)} 



(li 



IVc 



e=e 5 ; - 7=l,...,n 
BLOq 



u? LO °\an c 



e— ej ; j= l,...,n 

0, 



and for j = 1, . 



(36) 



A BLOj . BLOj 

= $ 



div^u^ 

BL0 3 



BLO, 



BLOj 



= -c^div | lx,(^)r i (iv 1 (6),o 

0. 



v f (x,(lf)lf 



if£eft 0j , 



rff + iV|X,(li e ), 



The variable ^ is opposite to £ / , i.e. to the first component of the vector Tj £ . So £ x 3 = Tj £ , where 
fj = IdTj and Id is the diagonal matrix with the diagonal elements —1,1. The constants c^ 3 , a, 3 are 
defined in such a way that the functions c\ 3 J^ 1 v(s)ds + d\ 3 and c\ 3 J^ 3 ^ Xl v(s)ds + d^ 3 are equal, i.e. 



(37) 



d 1 - c* 



v{s)ds + d e . 



Assume that every term in the sum in (|35[) is defined only in the branch of ilo 1 corresponding 

e— ej \j — l,...,n 

to e = ej, and it vanishes in 70. 

The solutions of these boundary layer problems decay exponentially at infinity and the constants c l 3 , 
c^ 3 , d\ 3 and d\ 3 are chosen from the conditions of existence of such solutions (see [5]). Let us define first 



c ( J from the condition of exponential decaying of u 



BLO, 



at infinity: 



^div£(fc-(gOrj(tfi(£),o) T )# 
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i.e. 



(38) - / i\M^)d£lC = / (tfgA de 2 S lfi , 

where the upper index 1 corresponds to the first component of the vector. 

Then we find c ; 3 , and as defined in ([37]). Then we determine the constants from the condition 
of the exponential decaying of pf LO ° at infinity. To this end, consider first problem (|3"5)) without the last 
term in equation (pJSjl -p i.e. 

-BLOo _ 



(39) 



£ [cf (^ a c (xM)) r, (mi®, of - v c (xMm) } 

I 

iv«cf (xM)rj (^(£2), of) , e e n O0 , 



div e sf LOn 



e—ej ;j=l, . . . ,n 

d 



£ 

e—ej ;j=l,...,n 



-BLOn 



^ |aoo = 0- 
Here the constants c ; 6j are just defined by (pIT|) and 



/ S 



e—ej ; j=l,...,n 



and satisfy the condition 
div^ (x,(^)r, (^Vife e ),0) T ) d£ 2 e = 



i.e. 

(40) 

Indeed, the choice of constants cf 



£ / ^(£1)^1 = 0. 

e=ej ;j=Q,l,...,nhj 

= — and cf J from (|38|) and condition (l29l) give relation ((40 



It is known that there exists the unique solution {u 



BLOo -BLOo 



Pi 



'} of this problem such that 



-BLOo 

h 



stabilizes to zero at infinity on every branch of Qq and pf LO ° stabilizes on every branch of Qq associated 
to ej, to its own constant pf L These constants are defined uniquely up to one common additive 

constant, which we fix here by a condition pf L °^ = 0. Then we define 



d l+l = -Pi 



(41) 



-BLOq 00 3 



BLO„ 



-BLOo 



BLOo _ f.BLO 
Pi — Pi 



e—ej ;j=X,.. . ,n 



on every branch of fig? associated with ej, i.e. pf = pf 



—BLOnooj 

*Pj Xj- 



Obviously, this pair {uf LOa ,pf LO °} satisfies (|55"j) . The boundary layer functions uf LOj and pf LOj 



j = 0, 1, . . . , n are not defined in the vicinity of O. Therefore we should change a little bit the formulas of 

, n. 

> ^ - f3 and 



u a and p a far from the nodes Oj, j = 0, 1, . 

Let r}j (x\ 3 ) be a smooth function defined on each segment ej , let it be one if x x 

let it be zero if x\ 3 — ^ < ^p. Let rj(x) = rjj{x\ 3 ) for each rectangle i3| and let 77 = 1 on each 7J. Set 

ryu) (x) = r)(x) on 7? and all half-rectangles having common points with 7? and extend it by zero on the 
remaining part of B|. Then we define u a and p a as 



(42) 



Z=0 



.i+2 



£ U f(^) Xe (x)+^^ 

e— ej ;j— 0,l,...,n i— 



BLOi 

h 



x - O, 



(43) 



p a = J2 £l+1 \ E Pf(* e,i )x 8 (aO + £Pi 

+£ £ i gf(^) Xs (z) + £g e ^ = 0) (1 - Xb (x)) *,(*). 

i=0 i=0 



77 w (z) 
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Let us mention that the last term (sum) in (|43|) corresponds to the boundary layer function pf LOj for 
I = -1. 



3.2 Error estimate 

In this section we estimate the error between the exact solution and the asymptotic one. Substituting the 
asymptotic expansions (|4"!S)) . (|43[) . into (|28p. we get the relations 



-div(i/(ir)£>u a ) + Vp a = f(x) 



(44) 



-v ? -xi(^)4+i)} + *> 

divu a = ip 
u a = g 



in B £ , 

in B £ , 
on <923 £ , 



where F e fc is the residual described in (fl~5j) . and -0 is defined by 



r o 



(45) 



-V x 77i(a;f ).u 



(BLO ) ( x ~ 



e 4 N „(BiOi 







in £ e n {sj* < /?} , 
m£ £ n{/3<af < 2J|ii}, 

inB e n |5^ <a .e i < | e; || ; 

in7f>« = 1, ...,n, 



ll^'ll(L 2 ((B e )) 2 = ^( ex P("r))' ll^'llfl' 1 (B e ) = O (exp(-^)) with a positive constant c and / ipd 



0, 



because / (u a ,n)ds — / (<7,n)ds = 0. The exponentially decaying residuals and ip appear from 

JdB e JdB' 

the truncation of the boundary layer terms by the function rj: it is different from 1 in the part of the 
domain B e where uf L i ( x ~f i ) , pf LOi ( z ~ 0i ) and their derivatives are exponentially small. We are going 
to prove the estimate 



l«-U lltfi(B«) 



= O 



We can not apply directly the a priori estimates because u a is not divergence free. 
Let us construct a function U a : B e — > M 2 satisfying the following properties: 



(46) 



div x U a = ip in B e , 

U a = in B e n {x\ < d e} \J^J = 1,.- 



Proposition 3.2. Problem (146|) /ias at /east one solution, satisfying 

\\U a \\ (H r m )*=0(exp(-c/e)) 

Proof. Due to (|4"6"jh we can consider the problem (|4*5]) as separate problem on each Bj. Denote 
by U ej the restriction of U a on S| , obviously J7 6j (x^ 3 , ) = for (x S S e such that < d e. 
For all (a; G S| such that x^ 3 > o?o£ introduce the new variable {y\ 3 , y^ 3 ) = ( r 1 - ^ j° £ i~f~); obviously 
iVi^l 1 ) S (0, 1) x (-5,5). Define a new function /x? : (0,1) x (-§,5) — >■ R 2 , by ^{vWvV) = 
^i J ((l e il ~~ do£)y\ 3 + doe, ey^ 3 ), jU^ ( ( | e j — d e)yl 3 + ^o£j sy^ ) ) ■ Obvious computations lead 



V |e, I -doe 

to the following problem for /x/ 



(47) 



diVjX" = "0((l e i I ~ d e)y e 1 3 + d e, ey e 2 3 ) in 
fit 3 =0 on 



(0,l)x(-i,i), 
fl((0,l)x(-|,i)), 
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(49) 



Applying the result of [3], Chap. Ill, p. 127 we get: there exist a solution of ITT)) such that 
( 48 ) M 3 ll (ffl((0 ,i)x(-i,i))) 2 - 0(exp(-c/e)). 

Expressing the norm of U ej with respect to the norm [i e 3 we obtain 

i-e- Il^ ej ll (ffl(f5 , n{ ^ <£i ~ 06})) 2 =0(exp(-c/e)). So, \\U a \\ (Him)2 = 0(exp(-c/s)). 
Define U a = u a — U a . Then (U a ,p a ) satisfies the following problem : 
-div{v(x)VU a ) + Vp a = f(x) + Aw{v{x)VU a ) + * 

£ { e fc+1 r, (f*) t xs(x) - £ fc (v a ,(^)< J + i - v^-(^')^i)} . in ^ 

e— ej ,j — l,... 7 n ^ ^ 

div« a = in B £ , 

u a = g on <9£> e , 

Theorem 3.1. Let (u a ,p a ) be the asymptotic solution given by (|42[) . (|43|) and (u £ ,p e ) the solution of 
([21]), £/ie following estimates hold: 

(50) f|| S --u«IW= (^ 

Applying a priori estimates and ([48)) we get 

l|C/ a -^|| H1(B .) =o(e fe+ - 

and so, 

||« -« e ||Hi (B .) = 0( E * 

The estimate for the pressure is obtained from the a priori estimates. These estimates justify the con- 
struction of the asymptotic expansion. 

Remark 3.1. The main result can be easily generalized in case when the length of (3j is different from 1. 

Remark 3.2. Formula (|38[) shows that only Cq 3 could be different from zero. The same analysis can be 
provided in the case of a multi-bundle structure, that is the union of a finite number of thin domains of 
B £ type (see section 4-5.2), and in this case the constants c( 3 should be determined from a system of 
linear algebraic equations (see (4-5.43), (4-5.44) in ISO])- 



4 Numerical experiments 

1. We consider the Stokes flow in a rectangular domain (0, 1) x (0,e) with e = 0.1: 



' -dxv(v(x 1 )Du) + Vp = Q in (0, 1) x (0,e) 



(51) 



divw = 

u{0,x 2 )=e 2 ^-(\ 



X2_ 

e 



tt(l,a 2 ) = e —\}-—) 



^ u(xi,0) = u(xi,e) 
here v{x\) = 2x\ + 2. From ([TO]) we get for j = 



in (0,1) x (0,e) 
on (0,e) 

on (0,e) 
on (0,1) 



r am %\ 

dxi \ v(x\) dxi J 

_ 2 1 dq x 2 



uq.i = £ 

"0.2 = 

Po = 



v(xi) 



_0qo_X2_ f _ X2\ 
l) dx\ £ V £ J 
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with the compatibility condition (|20p we have : 

/ 6(i-6)-"o.i(o,6)d6 = o 

Jo 

Since we assume that p(l, x%) = we get : 



q (xi) = -Xi(xi + 2) + 3 
2^ (1 _ ^ 



Uo,i = £~ 

UQ.2 = 

{ Po =0 



Solving numerically the problem (|5Tj) (by Comsol) we get the following results : for the first component 
of the velocity we have : 



Surface: Velocity magnitude (rWs) 




0.02 0.04 0.0S 



Figure 6: First component of the velocity Figure 7: The profile of the first component of the 

velocity for x\ = 0.5 
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and for the pressure 



Surface: Pressure (Pa) 




Figure 8: The pressure 



Figure 9: Pressure profile in x 2 = 0.05 



The error between the numerical result and the leading term of the asymptotic expansion is as follows 
We see that the error is of order of e et e 5 for the pressure and the velocity respectively. 
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Figure 10: Error for the pressure Figurc n . Error for the velocity 

2. Consider now the Stokes problem (J5TJ in a T-shape domain B e = (-1, 0] x (0, e) U (0, e) x (-0.45, 0.55) 



(52) 



-2 X 2 



t2 xi 

e 



t) 












■('■ 


Xi 





on (0,e) 
on (0,e) 
on (0, e) 



as inflow/outflow conditions, s = 0.1, u = anywhere else on the boundary and u(xi,X2) = 2 + 2(xi + 
1)/Lti(xi) + 2(cc2 + I)fi2(x2) + 2(^2 + l)/^3(x2). /ij are defined in such a way that v is equal to 2a;^ + 2 in 
each rectangle and constant near the nodes. 
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The functions \i\,\i2 and /X3 have the following shape : 



Figure 12: fii 



Figure 13: \x-i 



Figure 14: jj, 3 



So we get 



8/ 



Surface: Velocity magnitude (m/s) 



Figure 15: Thin domain £> e 



Figure 16: Velocity magnitude 



We do two crop sections, the first in £>f = (0, e) x (0, 0.55) and the second in £>f = (0, e) x (—0.45, 0), and 
the results arc as follows : 
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We see that the error is of order of e 5 as in the case of a rectangle and estimates. It confirms the 
theoretical prediction in section3. 
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